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We use the method of Balister, BoUobas and Walters [3] to give rigorous 99.9999% confidence 
intervals for the critical probabilities for site and bond percolation on the 11 Archimedean lattices. 
In our computer calculations, the emphasis is on simplicity and ease of verification, rather than 
obtaining the best possible results. Nevertheless, we obtain intervals of width at most 0.0005 in all 
cases. 



I. INTRODUCTION 

In this paper we study site and bond percolation on 
planar lattices, in particular the Archimedean lattices, in 
which all faces are regular polygons and all vertices are 
equivalent. The 11 Archimedean lattices are shown in 
Figure [U labelled with the notation of Griinbaum and 
Shephard [lOj: each lattice is represented by a sequence 
listing the numbers of sides of the faces meeting at a 
vertex, in cyclic order around that vertex. 




Square: (4*) Triangular: [2,^) Hexagonal: (6^) 




Kagome: (3,6,3,6) (3'',6) (3, 12^) 




(3,4,6,4) (4,6,12) (3^,42) 




(3^4,3,4) (4,8^) 



FIG. 1: The 11 Archimedean lattices 
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In site percolation on a graph A, each vertex, or site, 
of A is assigned a state, open or closed. In independent 
site percolation the states of the sites are independent, 
and each site is open with a certain probability p. The 
definitions for bond percolation are similar, except that 
is the edges, or bonds, of A that are assigned states. We 
shall write Pp for the corresponding probability measure, 
suppressing the dependence on A and on whether it is 
site or bond percolation that we consider. 

The basic question of percolation theory is 'when is 
there an infinite open cluster\ i.e., an infinite subgraph 
of A all of whose sites (for site percolation) or bonds (for 
bond percolation) are open. It is not hard to see that 
there is a certain 'critical probability' pc, such that for 
p < Pc there is never (i.e., with probability 0) an infi- 
nite open cluster, while for p > pc there always is. For 
this and other basic facts about percolation, see Grim- 
mett 0, or BoUobas and Riordan for example. When 
we wish to specify the lattice, and whether it is site or 
bond percolation that we are considering, then we write 
p^(A) or pb(A). 

The exact value of Pc is known in rather few cases: in 
1980, Kesten [ll| proved that p^{Z'^) = 1/2, where 
is the square lattice. Shortly afterwards |12l |. he proved 
that Pc = 1/2 also holds for site percolation on the tri- 
angular lattice T. Later, Wierman [36*1 used his 'sub- 
stitution' method to give rigorous proofs of the values 
pJ?(T) = 2 sin(7r/18) and pj?(ff ) = 1-2 sin(7r/18) for bond 
percolation on the triangular and hexagonal lattices re- 
spectively; these values had been obtained heuristically 
much earlier by Sykes and Essam [H, [s^l • There are two 
further values that may be easily derived from these: the 
(3, 6, 3, 6) or Kagome lattice K is the line graph of the 
hexagonal lattice, so pl{K) — p^[H) = 1 — 2sin(7r/18). 
Also, the (3, 12^) or extended Kagome lattice K^ is the 
line graph of the lattice obtained by subdividing each 
bond of H exactly once, so 

p^(i^+) = p^(i72) = ^/^R^ = (1 - 2sin(7r/18))'/'. 

These are the only critical probabilities known for 
Archimedean lattices. Indeed, it may be that the exact 
values of the other critical probabilities associated to the 
Archimedean lattices will never be known; they may sim- 
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ply be numbers that have no simpler descriptions than 
their definitions as critical probabilities. 

Given the dearth of exact results, it is not surpris- 
ing that much effort has been put into the estimation 
of critical probabilities. Almost all results in this area 
are of one of two types: (1) rigorous upper and/or lower 
bounds, and (2) heuristic estimates based on computer 
calculations. There are also a few heuristic derivations of 
conjectured exact results; we shall return to this briefly 
in Section |Vl Examples of (1) are the bounds obtained 
by Wierman [37, 38, 3^ (see also Parviainen and Wier- 
man [2^) using his substitution method. Even with con- 
siderable work on efficient algorithms and extensive com- 
puter calculations, it seems to be hard to obtain narrow 
intervals between rigorous upper and lower bounds: of 
the intervals listed in [2^, two have width a little under 
0.01, but many have width 0.1 or even 0.2. 

There are a small number of recent, more accurate, rig- 
orous results, including an intervals of width 0.00135 and 
0.0046 for site percolation on the (3, 12^) and Kagome 
lattices obtained by May and Wierman [2C|. 

Turning to (2), there are so many papers on this topic, 
going back to the 1960s, that it is impossible to attempt 
even a representative list. Let us mention a couple of ex- 
amples, however: forp^(Z^), Reynolds, Stanley and Klein 
[2a | reported the impressively accurate non-rigorous es- 
timate 0.5931 ± 0.0006 aheady in 1980. In 1986, Ziff and 
Sapoval [i^ gave the exceedingly accurate estimate of 
0.592745(2). In 1990 Yonezawa, Sakamoto and Hori [13 
gave estimates for several Archimedean lattices, with er- 
ror terms a little over 10^^. More recently, many very 
precise estimates have been given using quite sophisti- 
cated methods, for example by Suding and Ziff [s^l, New- 
man and Ziff [2^ and Parviainen [23] ■ It is very 
likely that these estimates are extremely accurate; er- 
rors of 'about ±3 X 10~^' are claimed in [s^], and even 
smaller errors in [2^ [H, However, these estimates 
come with no mathematical guarantees, and it is hard 
to be sure how accurate they really are. Although the- 
oretical error analysis is sometimes given (see Ziff and 
Newman [4^ . for example), this is certainly non-rigorous. 
Even assuming unproved results about the scaling limits 
of planar percolation models gives only the asymptotic 
behaviour of these errors; it does not allow us to say 
anything about the relationship between a finite number 
of data points and the true value of pc- Also, there is dis- 
agreement about even the asymptotic form of the errors 
in some cases (see Parviainen |24|). and there are several 
instances where earlier estimates have been contradicted 
by later ones. 

Surprisingly, it is possible to give a result intermedi- 
ate in nature between a rigorous bound and a heuristic 
estimate: one can prove that a certain (random) proce- 
dure generates a bound that is correct with probability 
at least 99.9999%, say; in other words, one can rigor- 
ously generate confidence intervals for critical probabil- 
ities. Such intervals are typically much narrower than 
the 100% bounds, although nothing like as narrow as the 



(claimed) uncertainties for heuristic estimates. Results of 
this kind were first proved by BoUobas and Stacey 0] in 
the context of oriented percolation, and then by Balister, 
BoUobas and Walters [3] in the context of (unoriented) 
continuum percolation. Here we use the method of the 
latter paper, which applies essentially 'as is' to percola- 
tion on 2-dimensional lattices, to obtain confidence inter- 
vals for the site and bond percolation critical probabili- 
ties for all 11 Archimedean lattices. Indeed such an in- 
terval for site percolation on the square lattice was given 
in [H, namely [0.5919,0.5935] (with a lower confidence 
of 99.99%). Here, with greater computational effort, we 
obtain narrower intervals. 



II. METHOD 
A. The mathematics 

The method of Balister, BoUobas and Walters Q 
is based on a simple application of the concept of 1- 
independent percolation (also known as 1- dependent per- 
colation). A bond percolation measure on a graph A, 
i.e., a measure on assignments of states to the bonds of 
A, is 1 -independent ii, whenever S and T are sets of bonds 
such that the graph distance from S to T is at least 1, the 
states of the bonds in S are independent from the states 
of the bonds in T. In other words, roughly speaking, the 
states of vertex-disjoint bonds are independent. Such 
measures arise naturally in percolation theory, in par- 
ticular in static renormalization arguments (see Section 
7.4 of Grimmett ^), and have been considered by many 
authors. Although the assumption of 1-independence is 
weaker than independence, it is strong enough to ensure 
percolation if the individual bonds are open with suffi- 
ciently high probability, as shown by the following lemma 
of Balister, BoUobas and Walters [3|. 

Lemma 1. Let ¥ be a 1-independent bond percolation 
measure on 1? in which each bond is open with probability 
at least pq = 0.8639. Then the probability that the origin 
lies in an infinite open cluster is positive. 

If the value of po is not important, then a weak form 
of Lemma [T] (with po replaced by some constant smaller 
than 1) is more or less immediate from first principles (see 
BoUobas and Riordan for example). It also follows 
from the very general results of Liggett, Schonmann and 
Stacey [TH comparing 1- (or k-) independent measures 
on general graphs with product measures. 

Starting from independent site or bond percolation on 
a lattice A C M^, there is a natural way to obtain a 1- 
independent bond percolation measure P on the square 
lattice 7?: given a 'scale parameter' s > 0, partition 
into disjoint s by s squares Sy, v ^ 1? . For each bond e 
of 1? let Rf. be the corresponding rectangle, so if e = uv 
then Re = Su U Sy. Let be some event that depends 
only on the states of the sites or bonds of A that lie 
within i?e, and take the bond e of to be open with 
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respect to P if and only if holds. Since the rectangles 
corresponding to vertex-disjoint bonds of are disjoint, 
this defines a 1-independent measure P. 

Suppose that the squares Sv and events Ee are chosen 
so that the following condition holds: 

whenever there is an infinite path 
V0V1V2 ■ ■ ■ such that Ee. holds for each = 
ViVi+i, there is an infinite open cluster in the ^^"^ 
original lattice 

and then p is chosen so that 

Pp{Ee) > 0.8639 for every bond e of Z^. (2) 

Then Lemma [1] implies that Pc < p. Indeed, we have 
already noted that P is 1-independent, so from ([2]) and 
Lemma [T] there is a positive probability that the origin 
is in an infinite P-open cluster. But then (as 7? is lo- 
cally finite) there is an infinite P-open path starting at 
the origin, so from ([T]) there is an infinite open cluster in 
the original percolation with positive probability. Hence 
P 1^ Pc, required. This type of argument, but in a 
qualitative form where the value of po is not important, 
provides one of the many easy ways of deducing Kesten's 
pJ?(Z^) = 1/2 result from a suitable 'sharp-threshold' re- 
sult; see Bollobas and Riordan 

Here we follow Balister, Bollobas and Walters @ in our 
choice for the event E^.- For w 6 Z^, let A^, denote the 
subgraph of A induced by the sites in Sy For each bond 
e = Ml! of Z^, let Ae by the subgraph of A induced by 
the sites in Re = SuU Sy. Let Eg be the event that each 
of A„ and A^, contains a unique largest open cluster with 
these clusters part of the same open cluster in A^. Here 
'largest' simply means containing the most sites. Note 
that Ee does depend only on the states of bonds or sites 
within Re, so we do obtain a 1-independent measure. 
Also, it is immediate that ([T]) is satisfied. 

To obtain an upper bound on pc, it remains only to 
find a pair {s,p) for which Q is satisfied. Note that 
it will suffice to check condition ^ for (usually) one or 
(occasionally) two bonds e of Z^: without changing the 
graph structure, we shall redraw all our lattices A so that 
the vertex set is a subset of Z^, and so that horizontal 
and vertical translations through some small integer C 
act as isomorphisms of A. For u = (a, 6), a,b £ we 
shall take 

Su — {{x, y) : sa < X < s{a + I), sb < y < s{b + 1)}, 

where s is a 'scale' parameter with C dividing s. Thus 
all squares are equivalent with respect to the lattice. 
Furthermore, for the lattices with an axis of symmetry, 
our new representation will have the line x — y as an 
axis of symmetry. This ensures that all rectangles Re are 
equivalent, so Pp(i?e) = Pp(-E^/) for all e, /. When there 
is no axis of symmetry, we have to consider one rectangle 
with each orientation. 

So far, we have only discussed upper bounds; this is 
because we can obtain lower bounds by bounding the 



critical probability for a related lattice from above. In- 
deed, given a planar lattice A, let A* be the usual planar 
dual of A, with one site for each face of A, and a bond e* 
for each bond e, joining the two sites of A* corresponding 
to the faces in which e lies. It is 'well known' that 

p;?(A)+p^(A*) = l. (3) 

Thus, to bound p^{A) fr om below we may bound pJ?(A*) 
from above. 

Note that while ^ is widely assumed to be true in 
great generality, it has only been proved under certain 
symmetry assumptions. Under very general conditions, 
the upper bound pj? (A) -fpi?(A*) < 1 follows immediately 
from Menshikov's Theorem |2l|. For the lower bound, 
one shows that it is not possible to have bond percola- 
tion in A at a parameter p and also bond percolation in 
A* at parameter 1 — p. For lattices (doubly periodic, lo- 
cally finite planar graphs) with rotational symmetry of 
some order fc > 4, there is a simple proof of the lower 
bound due to Zhang; see Lemma 11.12 of Grimmett 
where this argument is presented for Z^. Recently, Bol- 
lobas and Riordan 5] (see also Q) have pointed out that 
this argument can be easily adapted to lattices with ro- 
tational symmetry of any order k > 2. This is important 
here: all 11 Archimedean lattices have such rotational 
symmetry, but two, the lattices (3^,4,3,4) and (3'^, 4^), 
do not have rotational symmetry of higher order. Even 
more recently, Sheffield [29| has given a much more com- 
plicated argument that proves ([3]) for lattices without 
further symmetry assumptions. 

For site percolation, let A'* be the (in general non- 
planar) graph obtained from A by adding a bond between 
any two sites in the same face of A; we shall refer to A^ 
as the site dual of A. One has 

p^(A)+p^(AX) = l; 

the comments above about symmetry assumptions apply 
in this case also. 

B. The statistics 

For sufficiently small scale parameters s, it is possi- 
ble to find a p for which ^ holds by enumerating all 
possibilities for which sites/bonds in Re are open, and 
so writing Pp(i?e) as a polynomial in p. Needless to say, 
this is impractical and gives poor results in practice. The 
key idea of Balister, Bollobas and Walters Q is to use a 
statistical approach, obtaining confidence intervals with 
precisely calculated error probabilities instead of 100% 
upper bounds. Indeed, suppose that we have a random 
procedure A for generating a pair {sa,Pa)i and that one 
can prove that, with probability at least 99.9999%, the 
random pair produced is one for which Q holds. Then 
{—00, pa] is a (random, as always) one-sided 99.9999% 
confidence interval for pc (see below). Such a procedure 
A is very easy to define; again, we follow [3], with one 
small modification (and with different numbers). 
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Suppose that we have somehow 'guessed' values of the 
scale parameter s and percolation parameter p for which 
we expect that ¥p{Ee) is somewhat larger than 0.8639. 
We then generate N = 400 random simulations of the 
configuration within Re, and count the number m of 
them in which Ee holds. If Pp{Ee) — tt, then m has 
a binomial Bi(7V, tt) distribution with parameters N and 
TT. In particular, if tt < 0.8639, then 

P(m > 378) < P(Bi(400, 0.8639) > 378) 

= 1.1489- •■ X 10"^ < 10^V6- 

If our simulation does give m > 378, we can thus assert 
with very high confidence that tt > 0.8639, i.e., that ^ 
does hold, which, as noted above, implies Pc < p. 

This is the heart of the method of Balister, BoUobas 
and Walters Q (and also of the related method of Bol- 
lobas and Stacey Q): no matter how we arrive at our 
'guess' for s and p, provided we only perform one 'final' 
simulation, the simple inequality above shows that the 
probability that we assert an incorrect upper bound for 
Pc is at most 10~®/6. Note that we may be unlucky: 
if TO < 378, then we can assert only the trivial upper 
bound 1. In terms of the description above, our random 
procedure returns the guessed values {s,p) if m > 378, 
and the trivial pair (s, 1) otherwise. Note that we have 
no bound on the probability that we get 1 as an upper 
bound, but this does not matter for the argument that 
{—qo,pa\ is a 99.9999% confidence interval. Of course, 
to obtain useful results, we want to be reasonably sure 
that we will have m > 378, and this is where the careful 
choice of parameters comes in. 

Here we modify the method very slightly: the choice 
of the number 378 gives us individual error probabilities 
that are smaller than 10^^/6. Hence, we can perform up 
to three different runs with different parameters s and p 
(which may depend on the results of previous runs), and 
choose the best bound given by a successful run. It is 
still true that each run has at most a probability 10~®/6 
of producing an incorrect bound, so the probability that 
our final bound is incorrect is at most 10~^/2. Bearing in 
mind that the same applies to the lower bounds (realised 
as upper bounds on a dual critical probability), we still 
obtain 99.9999% confidence intervals. 

A small side note: since the lattice (3'*, 6) does not 
have an axis of symmetry, we ran our method in two 
directions, horizontally and vertically. This means that 
we wanted the individual error probability to be smaller 
than 10~®/12 which was satisfied by requiring at least 
379 successes in this case. 

There are two advantages to this method: it turns out 
to be slightly more efficient (based on heuristic calcu- 
lations). Bearing in mind that we can stop after one 
successful run, we can perform three runs each of which 
has a 90% chance of succeeding, say, more quickly than 
one run for the same p but a larger s that has 99.9% 
chance of succeeding. Secondly, if we are not very confi- 
dent of our guesses, after a failed first run we can choose 



more conservative parameters for the second and third 
runs (for example, keeping p fixed but increasing s), to 
be very sure of obtaining reasonable bounds in the end. 



C. Random number generation 

So far, we have assumed the availability of a suitable 
source of random numbers. In practice, one usually uses 
a pseudo-random number generator. This introduces a 
possible source of error: it could be that there is some 
pattern in the output of the generator that affects the 
results of the simulations. To minimize the likelihood of 
this we used the well known and well trusted MT19937 
"Mersenne Twister" generator developed by Matsumoto 
and Nishimura ^19*1, as updated in 2002. See their website 
l] for the source code and related literature. 

It would be very easy to modify our program to use 
other random number generators, or even a hardware 
generator. 

The selection of a random number generator for sim- 
ulations is often glossed over; here we emphasize this as 
it is important for our results: the only assumption in 
our results (that our procedure produces 99.9999% con- 
fidence intervals for pc) is that the random numbers used 
in the simulation may be treated as genuinely random. 



D. Choice of parameters 

In this subsection we outline the purely heuristic argu- 
ments we used to choose suitable parameters for running 
our final statistical tests. The correctness of the results 
does not depend on the correctness of these arguments. 
For this reason we allow ourselves to use consequences 
of the very widely believed but, except for one lattice, 
unproved conformal invariance conjecture. This conjec- 
ture of Aizenman and Langlands, Pouliot and Saint- 
Aubin [l3] states (among other things) that, for any 
planar lattice, after a suitable affine transformation, the 
limiting crossing probabilities for large regions are invari- 
ant under conformal mappings, and, more precisely, are 
given by Cardy's formula Q. For more details see Bol- 
lobas and Riordan for example. As shown by Smirnov 
and Werner [sH, building on work of Schramm [13] and 
Lawler, Schramm and Werner [3, fisl . \vf\ . this con- 
jecture, if true, enables the values of certain 'critical ex- 
ponents' to be calculated. Note that the conjecture has 
been proved, by Smirnov only for site percolation 
on the triangular lattice; for all other lattices it is still 
open. 

Fixing the percolation model under consideration, i.e., 
fixing the lattice A, and considering either bond or site 
percolation throughout, let pc be the appropriate critical 
probability, and set 

/(s,p)=Pp(i?e) 



5 



for one fixed bond e of Z^, noting tliat the definition of 
the squares and hence of the event depends on 
our scale parameter s. It is not hard to convince oneself 
that /(s,Pc) tends to some constant 0<a<lass— s-oo, 
although this does not obviously follow formally from the 
conformal invariance conjecture. 

Turning to the p-dependence of /(s,p), it is natural to 
guess that for fixed s, for p not too far from pc, the func- 
tion f{s,p) will roughly satisfy the differential equation 

±f(s,p)^Cis)fis,p)il-fis,p)), 
dp 

where C(s) is a constant depending on s (and on the 
lattice). For one thing, f{s,p) should decay exponen- 
tially, and approach 1 exponentially, as p moves away 
frompc- Also, by the Margulis-Russo formula, -^f{s,p) 
is exactly the expected number of sites/bonds that are 
pivotal for the event Ee, i.e., such that changing the state 
of this site/bond from closed to open or vice versa alters 
whether E^ holds. If a site (say) v is pivotal, then E^ 
must hold in the configuration with v open, and not hold 
with V closed, so it is reasonable to guess that for fixed 
s and a fixed site v, the probability that v is pivotal will 
be roughly proportional to Pp(i?e)(l — Pp{Ee)). 

Up to a constant factor, C(s) above is just times the 
probability that a 'typical' site (or bond) v is pivotal for 
Eg at p = pc- Roughly speaking, v is pivotal if and only 
if, when v is open, two open clusters of (linear) scale s 
are joined which, if v is closed, are separated by a path 
of linear scale s. Hence, the probability that v is pivotal 
should scale as s~"*, where a is the 'multi-chromatic 4- 
arm exponent'. Roughly speaking, a4 is defined as the 
scaling exponent of the probability that there are four 
disjoint paths Pi, P2, P3, Pi from v (or from 'near' v) to 
points at distance s from u, with Pi and P3 open, P2 
and P4 closed, with the endpoints of the Pi appearing in 
cyclic order. Assuming conformal invariance, from [31] 
we have 04 = (4^ — 1)/12 = 5/4, so we expect C(s) to 
scale as s"^^^/^ — s'^^^. 

Putting the above together, it is reasonable to expect 
the function f{s,p) to have approximately the form 



1 -|- exp(a — bs^/^{p — Pc)) ' 

for some constants a and 6 > that depend on the lattice. 
Our procedure for choosing the final parameters {s,p) to 
use is as follows: first, numerically estimate f{s,p) for a 
fixed small value sq of s (typically 72) and various values 
of p. Then fit the data with the function above to give 
a rough estimate of a and b. Then calculate values Pi/3, 
P2/3 oi p at which the formula predicts /(soiPi/s) = 1/3 
and f{so,p2/3) = 2/3. Next, run more extensive simula- 
tions to estimate f{so,Pi/3), ^^'^ use these two datapoints 
to calculate better estimates of a and b. The reason for 
this step is that we do not expect ^ to give a very ac- 
curate description of the shape of the curve f{s,p) with 
s fixed and p varying, particularly when p is far from pc, 



so we wish to extrapolate from consistently chosen points 
on this curve. 

Finally, we aim to choose a (large) s and a p close to 
Pc such that f{s,p) is approximately 0.957; this is be- 
cause P(Bi(400, 0.957) > 378) is close to 90%, so with 
these parameters we have a good enough chance of ob- 
taining a valid bound, bearing in mind that we can per- 
form three separate runs. Extrapolating ^ this far does 
not give very good results; experimentally, when ([4]) is 
about 0.945, or a little less, the true value of f{s,p) is 
large enough. Of course, the larger s is, the closer p can 
be taken to Pc- The exact values of s and p were chosen 
based on the amount of computer time available, and so 
that we obtained intervals of width at most 0.0005 in all 
cases. 



III. COMPUTATIONS AND RESULTS 

Although our final results are confidence intervals, we 
are aiming for rigorous confidence intervals, i.e., we must 
prove that, for each lattice, our procedure has probability 
at least 99.9999% of producing an interval containing the 
true value (assuming the random number generator we 
used is well behaved). The main practical consequence 
of this is that we must ensure that we evaluate Pp(£'e) 
for rectangles Re that fit together exactly in the manner 
required for the argument in Subsection III Al 

The first step is to transform each lattice so that trans- 
lations through some small constant C in the x- and y- 
directions act as isomorphisms. Such a representation of 
the lattice (4, 8^) with C = 8 is shown on the left of Fig- 
ure [U in this drawing, the vertex set consists of all points 
(x, y) G with x + y odd (for some of the lattices we 
use X + y even) . The white central portion of the figure 
shows a square region Su with scale parameter s = 8. 
It is this drawing of the lattice that we consider when 
defining Su, Re and Ee- 

Note that s must be a multiple of C, so that all squares 
Su induce isomorphic subgraphs of the lattice. For the 
lattices with mirror symmetry (all except for (3**, 6) and 
its bond- and site- duals), we choose a representation 
with the line a; = y or a; = —y as an axis of symmetry; a 
rectangle Re corresponding to a horizontal bond e may be 
mapped into a rectangle Rf corresponding to a vertical 
bond by a reflection in either of these lines, so this ensures 
that all rectangles Re induce isomorphic subgraphs of A; 
thus our program need only evaluate Fp{Ee) for one fixed 
(horizontal) bond of Z^. For the lattices without such 
symmetry, we run the same program on two drawings of 
the lattice, related by reflection in the line x ^ y; the 
horizontal rectangle considered for the second drawing 
corresponds to a vertical one in the first. 

Most of the representations we use are modifications of 
those shown in Figure 3 of Suding and Ziff ^3] , most of 
which have a horizontal axis of mirror symmetry. Since 
we want a diagonal axis here, we have rotated may of 
the representations by 45 degrees, obtaining a graph on 
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FIG. 2: The lattice (4, 8^) drawn with vertices a subset of Z^, 
in original form and in squashed form. 



points in 1? with x + y even. Our representations for all 
11 Archimedean lattices A and their planar duals A* are 
shown in the appendix. We have omitted any lattice for 
which the critical probability is known exactly. (In each 
case the site dual A^ is represented in the same way as 
A, but with additional bonds added to every face.) 

For those lattices represented with vertices in 1? with 
x-\-y even, it is computationally more efhcient to modify 
the representation to make it more compact, by mapping 
(x, y) to ([a;/2j , y), say. An example for the lattice (4, 8^) 
is shown on the right of Figure [51 Note that an sjr by 
s rectangle in the compact form corresponds to an s by 
s square in the original, where 1/r = 1/2 is the ratio 
by which we have squashed the lattice when compacti- 
fying it. In the program files, this is stored as the field 
RATIO for each lattice; this squashing is undone in the 
printJattice routines. 

The program perc.c, available from our website 
reads in the lattice, assigns states to the sites or bonds 
randomly, and then finds the largest open clusters in the 
left- and right- halves A„ and A^ of Ag, the subgraph of 
A induced by the sites in R^. (In fact, to avoid using too 
much memory, these two processes are done concurrently, 
see below for details) . Finally, it tests whether these open 
clusters are joined in Ae. The open clusters are found 
using a simple incremental algorithm that scans Su from 
the left and from the right. The method used to find 
the largest open cluster is (a simplified form of) that of 
Balister, Bollobas and Walters 0], and works as follows. 
We divide the square Su into strips which are narrow but 
are sufficiently wide that no edge jumps an entire strip 
(i.e., all the edges meeting the strip are entirely contained 
in the union of the strip and its two neighbouring strips) . 

We find the component structure of the open subgraph 
restricted to the two left most strips. Then we look at the 
next strip and find the new component structure formed. 
At each stage we have an equivalence relation on the 
vertices in a strip where two vertices are equivalent if 
they are in the same open cluster in the part of Su to 
the left of the current strip. We also keep track of the 
size of each of these clusters, and the size of the largest 
open cluster we have seen so far. When we get to the 



right hand edge of the square Su we know exactly which 
vertices (if any) in that strip are part of the largest open 
cluster of A„. 

We repeat the process on Su but working from right 
to left. Finally we add the edges between the right most 
strip of Su and the left most of Su and see whether the 
largest open clusters in each are joined. 

The important thing to note about this algorithm is 
that the storage required is proportional to the side 
length of Sui i.e. to s, not to the area of 5'„. 

IV. RESULTS 

For our percolation bounds see Table|Tl For full results, 
including numbers of successes, please see our website 0. 
Note that in the 400 simulations associated with each 
bound (or with each attempt to obtain a bound) we have 
always seeded the random number generator with 400 
consecutive seeds starting from 12345678. This means 
that the exact results of our simulations should be re- 
producible as a way of checking the program. Also, it 
shows that we have not performed many different runs 
and finally chosen seeds that work! 

The computations were performed running in the back- 
ground on around 70 (mostly fairly old) computers in 
the Department of Pure Mathematics and Mathemati- 
cal Statistics, University of Cambridge, over a period of 
around 2 weeks. This was made much easier by the fact 
that the department uses Linux rather than Windows! 
We are grateful to the computer officer, Andrew Aitchi- 
son, for technical assistance. 



V. CONCLUSIONS 

We have shown that it is practical to use the method 
of Balister, Bollobas and Walters 3] to obtain narrow 
confidence intervals for the critical probabilities for site 
and bond percolation on Archimedean lattices. Unlike 
the (presumably) much more precise estimates obtained 
by other methods, these intervals come with mathemat- 
ically guaranteed error bounds. The intervals are much 
narrower than those that can be 100% proved, and the 
error probabilities are very small; the running time does 
not increase much with a large decrease in the desired 
error probability, so a probability that is in practice zero 
(here 1 in a million for each lattice) may be achieved. 

We have tried to keep the computations relatively sim- 
ple; there is no point in using an algorithm that is proved 
correct if it is not possible to verify the computer pro- 
gram used. At the cost of more complicated programing, 
better results could be obtained in two ways. Firstly, the 
current program could be made to cache better and hence 
run faster by scanning the rectangle i?e in a more com- 
plicated manner: this 2s by s rectangle could be broken 
down into fc by fc squares small enough that the bound- 
ary of one square fits into the processors primary cache. 
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Lattice 


Site 


Width 


Bond 


Width 


Square 


[0.5925,0.5930] 


5 X 10"'' 


0.5 







Triangular 


0.5 





2sin(7r/18) 







Hexagonal 


[0.6968,0.6973] 


5 X 10"'' 


1 - 2sin(7r/18) 







Kagome 


1 - 2sin(7r/18) 





[0.52415,0.52465] 


5 X 10-^ 


-4 


(3, 12^) 


^1 - 2sin(7r/18) 





[0.7402,0.7407] 


5 X 10" 


-4 


(3,4,6,4) 


[0.6216,0.6221] 


5 X 10"'' 


[0.5246„0.5251] 


5 X 10" 


-4 


(3^42) 


[0.5500,0.5505] 


5 X 10"'' 


[0.4194,0.4199] 


5 X 10" 


-4 


(3^4, 3,4) 


[0.5506,0.55105] 


4.5 X 10"" 


[0.4139,0.4144] 


5 X 10" 


-4 


(3^6) 


[0.57925,0.57975] 


5 X 10"" 


[0.4341,0.4345] 


4 X 10" 


-4 


(4,6,12) 


[0.7476,0.7480] 


4 X 10"" 


[0.6935,0.6940] 


5 X 10" 


-4 


(4,8^) 


[0.7295,0.7300] 


5 X 10"" 


[0.6766,0.6770] 


4 X 10" 


-4 



TABLE L Rigorous 99.9999% confidence intervals for critical probabilities for site and bond percolation 



and these squares could then be processed column by col- 
umn. The overall storage requirement is approximately 
the same (one entire column must be stored) , but the fre- 
quency of cache misses is reduced by a factor of about k. 

A more significant improvement could be obtained by 
considering a different event : let i?e be the event that 
there is an open path crossing Re from left to right, and 
that there is an open path crossing the left-hand end 
square of Re from top to bottom. As noted by BoUobas 
and Riordan for example, this event still has the prop- 
erty ((!]). (Essentially this observation was used by Balis- 
ter, BoUobas and Walters [3] in obtaining a lower bound 
on the critical parameter for a certain continuum perco- 
lation model.) Also, the scaling behaviour of Pp(i?e) near 
Pc should be the same as for the event considered here. 
The gain is that whether or not Eg holds in a given config- 
uration can be tested faster, using an interface following 
algorithm of the type used by Ziff and Cummings [isj in 
1984, for example. Assuming conformal invariance, the 
expected length of the interface is s^""^ = s^^^, where 
— (3^ — 1)/12 is the multi-chromatic 3-arm exponent 
for SLEq. Note, however, that to use this algorithm in 
practice without running into memory/caching problems, 
one needs to generate the state of each site/bond from 
a pseudo-random function, rather than a pseudo-random 
number generator. 

All the exactly known critical probabilities associated 
to Archimedean lattices are roots of (simple) polynomial 
equations with integer coefficients. While it is easy to 
construct other lattices whose critical probabilities may 
be found in this form (see Ziff [i^, for example), it may 
well be that there are no such expressions for the re- 
maining Archimedean lattices, although some have been 
conjectured. In particular, Wu |4Qi] conjectured that for 
bond percolation on the Kagome lattice, pj? = 0.524429, 
a root of the equation — 6p^ + 12p^ — 6p^ — 3p^ 4-1 = 0. 
Tsallis [1^ conjectured the values pj? — 0.522372 and 
pj? = 0.739830 for bond percolation on the Kagome 
and (3, 12^) lattices, respectively. Tsallis's conjectures 
have been effectively ruled out some time ago by ex- 



perimental estimates (see [4l|, for example); they are 
rather far from the current best estimates of 0.5244053 
and 0.74042195. They have not yet been rigorously dis- 
proved, although the latest results of May and Wierman 
[2(j | come close. For both lattices, our results provide a 
rigorous '99.9999% disproof of Tsallis's values - they he 
outside our rigorous 99.9999% confidence intervals. 

Wu's conjectured value for p]!{K) seems to be much 
closer to the truth; it is well within the confidence in- 
terval we obtain. Nevertheless, it is still believed to be 
false; see Ziff and Sudin g [iq. for example. More re- 
cently, Scullard and Ziff [2^ have predicted certain val- 
ues for pj? for the Kagome and (3, 12^) lattices, using a 
heuristic version of the star-triangle transformation. Al- 
though they leave open the 'possibility' that one of these 
values might be exact, there seems no reason (to us, or, 
apparently, to them) to really believe this: the method 
is (as they admit) non-rigorous, and the value obtained 
in the same way for the Kagome lattice (given earlier 
by Hori and Kitahara without derivation) is outside the 
error bounds of existing experimental results. 



VI. APPENDIX 

Figure [3] shows the representations of the Archimedean 
lattices A that we used; the representations of the planar 
duals A* are shown in Figure ID The representations of 
the site duals A^ that we used are based on those of the 
original lattices A, with extra edges. The planar duals 
A* of lattices for which the critical probability for bond 
percolation is known exactly are omitted. 

Let us make some remarks about specific lattices. We 
wanted our lattices to have an axis of symmetry at 45 
degrees to the horizontal through a corner of each funda- 
mental region Su- in all the pictures except (3, 12^) this 
is the line down and to the right. In the bond dual of the 
(4, 6, 12) lattice the picture looks asymmetric but that is 
only due to our squashing to make it fit the square lat- 
tice. In other words if we reflect the graph about a line at 



Square: (4 ) Triangular: (3 ) Hexagonal: (6 ) Kagome: (3, 6, 3, 6) 

J/ 




(3*, 6) Horizontal (3'', 6) Vertical 




(3, 12^) 





(3^42) 



(3^4,3,4) 



(4,8^) 




(3,4,6,4) 



(4,6, 12) 



FIG. 3: The 11 Archimedean lattices on a square grid. 



45 degrees to the horizontal we get an isomorphic graph. 

One lattice, (3^, 6), does not have an axis of symmetry, 
so wc ran the program horizontally and vertically on this 
lattice. Both representations are shown and it is easy to 
see that one is the reflection of the other about a line 45 
degrees to the horizontal. 



For efSciency we tried to avoid having holes (vertices 
of the lattice not involved in the graph) in our represen- 
tations: it was not practical to avoid this for the Kagome 
and (3, 12^) lattices (and their site duals) and the bond 
dual of (4,6,12). 
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Kagome (3,6,3,6) (3'*,6) Horizontal 



(3*, 6) Vertical 




(3, 12^ 



(3^42) 



(3^,4,3,4) 




(4,8^ 



(3,4,6,4) 



(4,6,12) 



FIG. 4: The bond duals of 9 of the Archimedean lattices. 
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